The effective action and equations of motion of curved local and global vortices: Role 

of the field excitations. 



A. A. Kozhevmko\0 
Laboratory of Theoretical Physics, S. L. Sobolev Institute for Mathematics, 
and Novosibirsk State University, 630090, Novosibirsk, Russian Federation 
(Dated: October 6, 2008) 

The effective actions for both local and global curved vortices are derived, based on the deriva- 
tive expansion of the corresponding field theoretic actions of the nonrelativistic Abelian Higgs and 
Goldstone models. The role of excitations of the modulus and the phase of the scalar field and of 
the gauge field (the Bogolyubov- Anderson mode) emitted and reabsorbed by vortices is elucidated. 
In case of the local (gauge) magnetic vortex, they are necessary for cancellation of the long distance 
divergence when using the transverse form of the electric gauge field strength of the background 
field. In case of global vortex taking them into account results in the Greiter-Wilczek-Witten form 
of the effective action for the Goldstone mode. The expressions for transverse Magnus-like force 
and the vortex effective mass for both local and global vortices are found. The equations of motion 
of both type of vortices including the terms due to the field excitations are obtained and solved in 
cases of large and small contour displacements. 
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I. INTRODUCTION 



The active theoretical studies of the superfluid [J, [1| and superconducting vortices [3] take place over more than half 
a century Q. Important issues concerning the dynamics of both types of vortices include, in particular, the problems 
of the transverse Magnus- like force acting on vortex [1, 0] , and the evaluation of the vortex effective mass @, H, 0] . 
The above incomplete list of references shows that different models are used in describing various aspects of vortex 
dynamics. This is due to complexity of the object under study put in real experimental conditions which include 
dissipation, finite-temperature effects etc. The common feature of both types of vortices in condensed matter physics 
is that they appear as topological defects in the models with spontaneously broken U(l) symmetry. 

The purpose of the present paper is to study the dynamics of the curved line defects in the framework of a uniform 
approach based on the field theoretic models and to scrutinize the role of the excitations of the gauge and scalar 
fields exchanged between the different segments of the vortex. The local (gauge) vortex equations of motion are 
obtained from the effective action of the Abelian Higgs model (AHM) at zero temperature T = 0. The same task is 
performed for the global vortex by taking the limit of the vanishing gauge coupling constant. The earlier attempt to 
consider the role of the field excitations was undertaken in Ref. However, the background contribution which is 
absolutely essential for the static situation was not taken into account in Furthermore, the treatment in Ref. Q 
was restricted to considering only quasi-two-dimensional situation (straight vortices), and no attempt was undertaken 
to derive the vortex equations of motion by means of the variation of effective action. In the present work we fill these 
gaps and treat the background fields and the excitations. We work in truly three-dimensional situation and allow for 
the curvature of the vortex contour by deriving expressions for the effective actions which systematically use the so 
called derivative expansion. The expressions for the transverse Magnus-like force and the effective masses of the local 
and global vortices are obtained. 



II. FORMULATION OF THE MODEL AND BASIC NOTATIONS. 



We will work in the London limit where the London penetration depth is much greater than the coherence length 
^. [See Eq. (|3.2p and (|4.1ip below for the definition of these quantities.] In this limit the gauge field configuration of 



2 



the magnetic vortex is generated by the singular phase Xs of the order parameter tfj jlOf . The notion singular means 
that 

[Vx V]xs = 27r^na f daaX'j(-^\x ~ Xa), (2.1) 

a '' 

where sum goes over all contours given by the vector Xa = Xa{t,aa), Ua is the number of flux quanta trapped by 
the vortex. Hereafter, prime means the differentiation with respect to the corresponding contour parameter CTo, and 
[a X b] stands for the vector product of two vectors a and b. Since the length of the contour is / = J^°^ \X'^\d(Ta, the 
natural choice of aa is that it gives the length of the contour segment starting with initial value aai and ending with 
CTa- Hence, the gauge condition 

= 1 (2.2) 

can be applied. In the rest of the paper the limits of integration over contour parameter are omitted for brevity. The 
action of the Abelian Higgs model is 

d^xi^iE'' - H^) - |(|V^|2 - no)2 + ]p,*(iWt -q^ + qao)^ + c.c.]- 



2m 



-inV - -A + -a] tp 

c c 



poifi } . (2.3) 



Here, E = —dtA/c— 'Vip and H = 'V x A are the electric and magnetic field strengths, with (p standing for the time 
component of the gauge four- vector potential = ((^, A); q = 2e, m = 2rne are the charge, mass, of the scalar field 
in terms of the electronic ones, c is the velocity of light. Further, uq and po stand for the density of the scalar field 
condensate (the number density of the Cooper pairs) and for the homogenous positive charge density introduced to 
provide the net neutrality of the system: 

Pq + qno ^ 0. (2.4) 

The coupling constant g will be related below to the sound velocity c^. The four- vector = — — dfj,Xs represents the 
four-gradient of the singular phase. The field strength tensor corresponding to is /^^ — d^a^, — d^a^ = {Eg, Hg), 
where the mixed Fourier transforms of the singular parts of the electric and magnetic field strengths are 

Eg{t,k) = -^^nay"da,[X, xX;]e-''=^", 

Hg{t,k) = <foJ2na f daaX'^e-'"''^ . (2.5) 

a •' 

Hereafter dot means the differentiation with respect to time, and 

2'Khc 



$0 = 



is the magnetic flux quantum. Both the full Fourier representation in A: = (w, k) or the mixed one in (t, k) turns out 
to be very suitable in what follows. The full representation is defined as usual: 

bk = b{Lu,k) ^ J d^xe"^'-''''''b{t,x), 

b{t,x) = I (2.6) 

The mixed one is defined accordingly. One can find the necessary expressions for the singular part of the vector 
potential ak{t) = a{t, k) from the equation i[k x a^] = Hs{t, k): 

ak(t) ^ i^Y.'^a j daa[k X X'Je'^"''^. (2.7) 
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This vector is transverse, ka^ = 0. The mixed Fourier component of the singular part of the Coulomb potential oq 
can be found from the expression Esit, k) = —a^/c — ikaQp.: 

aofcW = -*^I]^a / d(Jak[Xa X X'Je-^'^''^ . (2.8) 

a •' 

Note that both E{t, k) and aofc look like the local Lorenz transforms of H{t, k) and a^, respectively. The end point 
contribution vanishes identically for the closed contour, while for the open one the end point term is proportional to 
the combination Xj ■'V\x — Xf\^^ — Xi • V|a; — which can be made vanishing either by pinning the end points 

Xf^i or pushing them to the spatial infinity. In what follows we will always assume the vanishing of the analogous 
end point contributions. 



III. THE CONTRIBUTION OF BACKGROUND FIELD . 



Let us rewrite the action Eq. (|2.3p . first, in terms of the modulus and the phase of the scalar field ijj = n^^'^e''^. 
Second, expand the field configuration in terms of the background fields and the fiuctuations, Lp — > Lp+5Lp, A — s- A-\-5A, 
n — > no + 5n, x ~^ ^X- Note that the space-time derivatives of the singular background phase is already included in 



Eq. (|2.3p as the four- vector a^. The modulus of the scalar field is taken to be no as it should in the London limit. 

as the sum of the background action Sbg and the action of the fluctuations S{. 
The background action is 



One gets the total action Eq, 
The latter is discussed in Sec. HVl 



S'bg - 



1 



{A - of + noqao 



(3.1) 



The term proportional to (p drops drops in view of Eq. (|2.4p . Hereafter 



Xl = 



(3.2) 



is the the London penetration depth. The condition of the vanishing of terms linear in SA and Sip gives the gauge 
field (Maxwell) equations for background gauge fields 



[V X H] 



1 



1 



1 



dtE, 



(3.3) 



= 4n{nnq + po) = 0, (3.4) 

while the vanishing of terms linear in variation of the phase Sx gives the equation 

V(A-a)=0. (3.5) 

The latter equation is satisfied automatically in the Coulomb gauge. Since the situation is non-relativistic, we will 
neglect the retardation effects in what follows. Hence, the displacement current in Eq. p.3p will be neglected, too. 
One should not vary over density because we assume homogeneous density n — uq everywhere except the vortex core. 
The latter is ignored in the London limit. Using the relation — i[k x A^] one can find the solution of Eq. (|3.3p 
in the mixed {t, k) Fourier representation: 



1/Ai ^ _ z$o/Ai 



k^ + i/xi " k^k^ + 1/XI) 



a J daa[k X X'Je-'"''^ . (3.6) 



Making the inverse Fourier transform to the coordinate space one can see that the vector potential A is decomposed 
into the long range and the short range terms, while the combination 



«fe = -TIT— ^ E / ^'^-['^ ^ X;]e^^'=^- (3.7) 
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which appears in Eq. (|3.ip contains only the short range exponentially damped contribution. The expression for the 
Fourier amplitude of the strength of the background magnetic field is 



\2 

The Fourier amplitude of the strength of the electric field is 



a J daaX'^e-'"^-. (3.8) 



l/(cAi) . 

The scalar potential (pk = 0, as can be found from the above equation and the Gauss law Eq. (|3.4p . The resulting 
expression for Ek looks like 



2 



fe2 + l/A' 

Both Hk and are transverse. As for it is evident from Eq. (|3.9p . In the case of Jf^ the transverse character 
is proved by noting that kH^ oc J dakX' e^^'^-^ = i f da-^e^^''-^ — for closed contours. For the open ones the 
transverse character is provided by pushing the end points to the spatial infinity. 
With the help of the relation 

(^I^WI' (3-10) 

one can obtain the action of background fields in terms of the vortex contour variable Xa = Xa{t,(Ja)- From now 
on we will restrict ourselves by the case of the single contour with the unit flux quantum. The notations adopted in 
what follows are Xi,2 = X{<ti^2), X12 = Xi — X2, where the points ai^2 belong to the same contour. Taking into 
account Eqs. ((231), dSH), (0121), and 0JU\i one obtains 




2 

noQ J d'^xaQ. (3-11) 

In the braces, the terms with the time derivatives originate from E'^, the term cx A| does from the kinetic energy of 
the Higgs field cx (A — a)^, while the remaining one is due to H^. Notice that the contribution of the term with oo 
in Eq. p. lip , at first sight, seems to be irrelevant because 

d^xaa = -— / dalX x X'] ( [ d^xV- ^ -] = 0. 

AttcJ \J \x-X\J 

However, its variation over X{a,t) is not zero and, as is found in Sec. IVI[ results in the Magnus- like transverse force. 

Let us obtain the leading and next-to-leading terms in the expansion over small parameter K^Xj^, where k = \X"\ 
is the curvature of the contour. The dominant contributions come from the nearby points of the contour (T2 = (Ti + z, 
z <C (Ji, hence the derivative expansion 

X2 = X,+ zX[ + + ^X'l' + • • • , 

2 

X21 = zX[ + -jX'i + y-^i" H 1 

X', = X[+zX'l+'^X'l' + ---, 

z2 



= |z|(l-_«2 + ...] (3.12) 
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can be used. When obtaining the bottom hne in Eq. ((3A2)) the relations X' ■ X'" = -X"^ = -k^ and X' ■ X" = 
should be taken into account. They follow from the Frcnct-Serre equations 

X" = K,n, 
n' = -kX' + rb, 

y = -TTi, (3.13) 

where t stands for the torsion of the contour, and n, b are the vectors of normal and bi-normal, respectively. The 
vectors (n, b, X') comprise the right triple of the unit orthogonal vectors, so that X' = [n x b] (and similar relations 
obtained by the cyclic permutation). Note also that X2 = 1 + 0{z'^) if X'j^ — 1, hence the approximate expansions 
over z do not break the gauge condition X'^ = 1 within the adopted accuracy. The background action is represented 
as ^bg — sl^J + AShg, where 

$2 r ^ d'k f l/Ai 



/ ( k^ + l/Ai ) / ^'^i^'^^e-^^- {[Xi X X[] ■ [X, X X',]/c' {X[X',){1 + \lk')} + 



noq J d'^^xao- (3.14) 
It is finite at large distances while 



(27r)3fc2 \^k2 + l/A2 



j daida2e-'^^''-{k[Xi x X[]){k[X2 x X'^]) (3.15) 



diverges at large distances in view of the factor 1/fc^. As it will be shown below, the similar divergent contribution 
with the opposite sign arises from the term originating form the integrated out fluctuations. Hence, we postpone the 
evaluation of the terms cx l/fc^ to Sec. fVl 

Performing the integration over momentum k in Eq. p.l4p one gets 



The integrand in Eq. (|3.16p is nonlocal in the contour parameter a. This is due to the distribution of the gauge 
vortex profile over coordinate space. However, the contributions of the remote points are suppressed exponentially as 
g-|x((T2)-x((T2)|/Ai, ^ gj^f^ (.g^j^ dcrlvc the approximate local expression. Indeed, using /dui / da2 ~ / rffi f^^dz 
and replacing ai a one can integrate over z to obtain 

/ \ 2 /- 



nog j d'^xao- (3.17) 

where C = 0.577215.... is the Euler constant. Notice that the term with c>c InA;/^ arises due to the usual short 
distance cutoff when the lower integration limit over z is replaced by the coherence length ^ : 



oo 



dz-^^ >2 ^e-"/^^»2 In^-C . (3.18) 

\A Ji z \ i J 

Such a replacement is justifiable because the modulus of the scalar field 1-01, in fact, vanishes at the vortex core rj^ < ^. 
Thus, the only contribution to the action enhanced as InA^/^ is that of the kinetic energy of the scalar field given 
by the term cx k^ in Eq. p. lip . It is recognized as the contribution of the minus energy of the gauge vortex with the 
energy per unit length ey given, in the leading logarithmic approximation and by neglecting higher derivative term 
(X K^, by the expression [lO| 

^° ^'in^. (3.19) 



47rALy e 

One should be cautious when dealing with the approximate local form Eq. (|3.17p [together with the correction given 
by Eq. (j5.10p ]. Such local expressions are useful only when identifying various contributions, as, for example, the 



energy per unit length above, or in deriving the expression for the effective mass of the vortex (see Sec. IVIl IVIIi and 
Ref. [9| ) . They cannot be used for obtaining equations of the vortex motion because of the above mentioned finite 
size vortex gauge field distribution cx e^\^i'^2)~^i'^2)\/XL ^ Indeed, the variation of the approximate local expression 
oc / dtdaX'^ in Eq. (|3.17p gives the result which is by the factor two greater than the correct expression obtained 
upon varying the corresponding nonlocal contribution (x / daida2X[X2e~^-^^^^^'^'^ in Eq. (|3.16p . The same is true 
for the dynamical terms containing the velocity X. See Sec. IVIl for more detail. To summarize, only expressions 
Eq. (|3.16p [and (|5.8[) in Sec. IVI| should be used to obtain the vortex equations of motion. 



IV. INTEGRATING OUT THE EXCITATIONS OF THE SCALAR AND GAUGE FIELDS. 



The next step is to integrate out the propagating fluctuations of the modulus of the order parameter Sn, the phase 
of this parameter Sx, the Coulomb potential 6ip, and the vector potential 6A. This is necessary for obtaining the 
correction to the action of the background fields. To this end let us write the action Sf for fluctuations representing 
them for short as the row f'^ = {Sn, (5%, Sip) and SA. One obtains Sf — 



(0) 



where 



:.(o) _ 



1 



1 



and 



S, 



(1) 



Stt 



qh 



-J' Mf - —6A -^di + V' ^ -^\5A + 5' 



{A -a) 



The Coulomb gauge V(5A = is chosen. The matrix M is 



4mno 




(4.1) 



(4.2) 



(4.3) 



Formally, S^^"^ in Eq. (|4.2p is quadratic in fluctuations and should be attributed to the free action. However, because 
of highly inhomogeneous profile of the background vector potential A — a, the path integrations with S[ cannot be 
performed in the closed form, so we prefer to treat S^^^ as the perturbation leaving 5"^°^ in Eq. (|4.ip as the action 
of free fluctuations. The contribution of the term q^{A — a)^/2mc^ can be neglected. In fact, it would enter the 
effective action cither quadratically giving the terms oc J J J J daida2da3da4X[X2X^X^, or as the interference with 
quo giving the terms (x J J J daida2da^XiX[X2X'.^. Both type of the terms are attributed to higher derivatives and 

go beyond the scope of the present paper. One can show also that keeping sf'^ results in the higher derivative terms 
in the effective action. See the end of the present section. 

As is known, one should find the generating functional in order to integrate out the quantum excitations 
In the present case it looks as 



I 



D[5n]D[5x]D[5ip\D[5A] expi-fi-f] 



(4.4) 



where the action Si = S\ + ^ d'^x{jj f +jASA) is the functional of currents = (jn, jx' Jv)' ^^"^ i^- The gauge 
fixing and ghost terms are omitted because they are irrelevant in the context of the paper. As in the preceding section, 
it is suitable to work in the momentum space k = (uj,k). Then the integration measure for the spacetime function 
f{x) is expressed in terms of its Fourier amplitudes as 



Dif] - n 



dRefkdlmfk 
2tt 



The path integral Eq. (|4.4p is gaussian, so using the inverse of M in the momentum space, 
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one finds the generating functional Z[j] = exp j-AS{^ where the correction due to excitations is 

~ 1 f d^k 

ASi = i?ilndetM + - 




^2 172-^ ] + —[Un + qao)Uj^)k- c.c.y 



'hk 



47r|(jA)s 



1 



(4.6) 



Here, 



+?^'c^fc^ (4.7) 

2m 



is the square of the Bogolyubov spectrum [l^] looking at small wave numbers as the spectrum of sound waves with 
the sound velocity 



^. (4.8) 

TO 

As is known, the regularization prescription — iO guarantees that the causal (Fcynman) Green functions will result 
after taking the variational derivatives of Z[j] over j. The correlation functions in the momentum space which are 
just the Fourier transforms of the Green functions are 

{f^H^ft ■■■) = Z-'n (-^^^Tl^)i f-2*^-i) ) • • • Z[i\\j=o- (4.9) 



5j 



The divergent constant In det M drops from all expressions for the correlators Eq. (|4.9p . 

The meaning of the poles in oj in Eq. (|4.5p is the following. Let us consider the freely propagating fluctuations 
decoupled from their source represented by the term in square brackets of Eq. (j4.ip . Their equations of motion 
obtained upon the condition of the vanishing variational derivatives look as 

= hdtSx + gSn + qSip — ^^Sn, 

= dtSn+^VHx, 
m 

= V^Sip + 4Trq6n. (4.10) 
Applying to the first line in Eq. (|4.10p and using the second and third lines of the same equation one finds 

The propagating plane wave solution has the dispersion law with the gap 



, 4(fc) 



which is just the pole position of the matrix M^^. The same dispersion law is obtained for the phase fluctuations Sx- 
The above dispersion law in the long wave limit acquires the form uj'^ u!p + c^k^ = oj^ + ^Vpk'^ , where usp, vp stand 
for the plasma frequency and Fermi velocity, respectively. This is the dispersion law for the Bogolyubov- Anderson 
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mode in the charged Fermi system [Ij, [Tj] . The sound velocity Cg is related to vp through the expression for the sound 

velocity c^s — m (^) ^•^ taking into account the expression for chemical potential of the weakly interacting 

Bose gas /i — gno/m [l^ Recall that AHM mimics the Bose liquid of the Cooper pairs whose parameters are 
related to its electronic counterparts by the relations /i = 2(1^, rig — m — 2m^^ q = 2e. In the case of neutral 

superfluid where the gauge coupling constant is switched off, q — > 0, the spectrum is gapless, u)^ = e%/h^ On 
the other hand, taking the static limit and neglecting the gauge fields one finds the equation 

V On T — on = 0, 

which means that the space fluctuations of the modulus of the scalar field damp at the coherence (or healing) length 

' (4.11) 



2{mnogy/'^ 

Taking into account Eq. (|4.8p one obtains the relation 

e = ^ (4.12) 

Zmcs 

which relates the sound velocity Cg with the coherence length ^. We dwell upon in this paragraph on the well known 
issues in order to show that the character of excitations and their dispersion laws can be established in the framework 
of effective lagrangians without references to the underlying microscopic picture. 

The contribution of the term Eq. (|4.2p can be evaluated with the help of Eq. (|4.6p . One finds the correction to the 
effective action due to the above term: 



Ac(i) "0 [ d'^kid'^ka . . , ,2 /eK^i) , 2 -nA 

^^ Am^n^ J (27r)8 (^ " "ifci+fc^^ifeifeaaofciQofc, ^2 + A| ^ ~ j ^ 



— — iO 

L 



Using Eq. p.8|) and (|3.7p one can see that because of the factor uji in the integrand the above expression starts with 
the terms of the type ^ diJidij2da3X'^X2X^X'^ etc. resulting in the third order time derivative terms in the vortex 
equations of motion. Such higher derivative terms can be neglected for sufficiently smooth evolution of the contours. 
By this reason the term Eq. (|4.2p can be neglected. 



V. DERIVATIVE EXPANSION AND THE EFFECTIVE ACTION FOR THE LOCAL VORTEX. 



The correction due to fluctuations is obtained from Eq. (|4.6p by setting jr, 
term Indet M: 



0,jx — — ^ii'i neglecting the 



AS'f 



d'^k 
(2^ 



-1 2 
noq 



\aQk 



2 1,2 



(5.1) 



where oofc = dte'^'^^ a^kit) , with aQk{t) given by Eq. (j2.8p . is the four-dimensional (w, k) Fourier transform of the 
time component of the singular part of the gauge field potential a^. The explicit expression will be obtained in the 
case of relatively slow dynamics of the vortex contour. In Sec. IIIH the derivative expansion based on the smooth 
contour shape is used to obtain the explicit expression for the background part of the effective action. In the present 
section, the analogous derivative expansion based on the slow motion of the contour is the key point. Let us expose 
it in detail taking the typical integral: 



duj 
2^ 



\fk 



iO 



dtdr 



duj Ju{t)fl{t + T)e'' 



dt 



27r L 



1 

2n 



iO 



= / dtfkit) 



00 

E 

(=0 



l\ dt^ 



1 



(5.2) 
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where dots stand for the terms with the higher time derivatives. In the case of our interest, a typical fk contains 
X, hence d^f^/dP, and, consequently, the equations of motion would contain the higher derivative term d'^X/dt^. 
In order to have the equations of motion with the time derivative not higher that two one should keep only the first 
term in the expansion Eq. (|5.2p . Hence, in the lowest order in the number of the derivatives over time, one finds from 
Eq. (|5.ip the correction to the effective action due to the fluctuations: 



ASf = 



8nXl 



dtdaida2 



{k[Xi X X[]){k[X2 X X'^]) 



fc2 



(5.3) 



One can see that due to the factor l/fe^ this expression is logarithmically divergent at large a2 — cri. It is the same 
large distance divergence as is observed in the background action Eq. (|3.1ip , (|3.f 5p , but with the opposite sign. Let us 
examine the divergent and finite contributions in Eq. (|3.I5p and (|5.3p . Denoting the their sum as AS'off = A^bg + AS'f 
and taking into account the relation Ara^c^g/h^ = I/^^, let us represent AS'cff in the form suitable for integration over 
momentum: 



A^eff 

/(cti, 0-2) 



dtda2da2[Xi x X[]i[X2 x X'^],"^ 211"^ 2ijl[(^i, (^2), 



Ai 



d^k 

(27r)3 



ikX 



(27r)3 
1 



1 

fc2 



1 



1/A2 



fc2 + 1/^2 



AhX2 



fc2 fe4_^fc7^2_^ 1/^2^2 



(5.4) 



where V2ii = d/ 8X211- The above equation displays clearly the cancellation of the divergent contributions arising at 
|fc| 0. Notice the appearance of additional length scale 



As = Al • - < Ai, 

c 



(5.5) 



see Eq. p.2p . (|4.8p . and (I4.12p . This scale is absent in the static case, but arises in the dynamical situation. After 
convincing oneself that the divergences cancel one can write the resulting expression for the quantity I{ai,a2) in 
Eq. (15^ : 



/(a'i,CT2) 



4^|X2i| 



1 



1^21 1 

1/2 



1 



87r|X; 



211 



v/1 - 4eVA^ 



-|X2i|(l + Vl-4?VA?}'^'/?\/2 



1 - 



VI - 4eVA? 



(5.6) 



at ^ < As/2, and 



/(0'1,CT2) 



^-\X2i\/\l 

4^1X21 1 



1 



1^21 1 

2Al / 4^1X21 1 



-|X2i|Vl/2«A, + l/4e2 



cos(|X2i|Vl/2eA. + l/4e^ 



sin 1X21 1 v/l/2eAs + 1/4^2 



(5.7) 



at ^ > As/2. Both above expressions are interrelated by the analytical continuation. The scale As divides the London 
part of the parameter space ^ <^ Xl into two pieces. As < £, Xl and f < As <C Ai,. Let us consider the latter one 
and, in addition, choose a special case ^ <C As/2 <C Al, leaving the case f > As/2 for a future work. Expanding square 
root in Eq. (|5.6p . neglecting the terms exponentially small at IX21I > ^, and applying to the resulting expression 
V2iiV2ij one finds 



AS-eff = 



[ 

327r2c2 J 



dtdaida2[Xi x X[],[X2 x X^ { - 



5i' 



\X2l\ 



1 



\X2l\ 
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-|X2i|/A, _ ^ 



1^21 1 



AX2i\/\l 



l^2l| 



A, 



l^2l| ^21 



A. 



2Ai|X2i| 

-|^2l|/Ae„ 



^-\X2i\/\l_ 



^21 

2Ai 



1 



1^21 1 
Ar 



HX2i|/At 



(5.8) 



Note that 5ij5^{X2i) arising due to V2ijV2ij (I/IX21I) is multiplied by the factor vanishing at X21 = and hence 
does not contribute to the effective action. The derivative expansion adopted in the present paper, results in setting 
(T2 = iJi + z and expanding all quantities depending on (T2 in series over z, see Eq. p.l2p . In view of the exponential 
damping of the integrands at large |-X"2i| one can replace in Eq. (jS.Sp approximately J dui J d<T2 = / d<Ti dz. The 
contributions proportional to 6ij possess the logarithmic divergence at short distances. This is an artifact due to the 
ignorance of the vortex core IX21I < One can handle this divergence in the same manner as in Eq. (|3.18p . Taking 
into account another regularized expression 



dz f z 



2 [e 



C-lnA/C 
A2 



(5.9) 



where C = 0.577215.... is the Euler constant, one can see that the terms oc ln(AL/C) cancel, the term oc ln(As/^) 
survives. Gathering all the integral together and replacing cti — > one obtains the correction to the effective action due 
to the combined contributions of the excitations integrated out, and of the part of the background action Ec^. (13. 15^ : 
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[X X X'] 
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{X[X' X X"]f 



(5.10) 



The higher derivative terms do not have the logarithmic enhancement factors InXj^/^ nor InAs/^. Hence, we neglect 
them to obtain 



\47rAi 



dtda IJ^[X X Xf\n^ - X'^ In^"- ^ ' > -^^^ 



i2cr 



+ TiQq I d xaQ 



(5.11) 



The contribution of the short range part of the electric field has the same form [X x X']^, see Eq. (|3.17p . But it is 
not enhanced as logarithm and is suppressed by the factor c^/c^. Hence, it can be neglected. Since [X x X']^ = X^ 
in the gauge X • X' = 0, X'^ = 1, one can see that the first term in the braces of Eq. (|5.1ip can be interpreted as 
the kinetic energy of the vortex motion with the effective mass per length L 



TOoff 
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A. 



(5.12) 



The first equality in Eq. (|5.12p coincides with the expression obtained earlier for the straight vortex [9|. The last 
equality shows that the effective mass of the gauge vortex equals to the mass of superconducting electrons expelled 
from the region with the transverse dimension ^, multiplied by the dynamical enhancement factor. 



VI. THE LOCAL VORTEX EQUATIONS OF MOTION. 



First of all, let us show that the variation of the term n^q J d'^xao in Eq. p.l6p gives nonzero result despite the fact 
that the contribution of the above term to the effective action vanishes. Using the Fourier amplitude of the singular 
part of the time component of the vector potential qq [see Eq. ()2.8p ] one obtains 

S J d^xao = J dtda{6X[X x X']). (6.1) 

Here, a number of integrations by parts over both t and a should be performed, and the known vector relation 
[Ax [B X C]] = B{AC) — C{AB) should be taken into account. The equations of motion of the local (gauge) vortex 
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in the approximation of large logarithms can be obtained upon varying Eq. (j3.16p and (|5.8p over the contour variable 
X and with the neglect of the higher derivative terms. To do so, one should first vary over X including the factor 
g-|X2i|/At ii^Qj^^ second, use the expansions Eq. (|3.12p keeping the terms with the lowest non-vanishing order in the 
expansion variable z = (T2 ~ en : 



6 / daida2{X[X'2] 



-\X2i\/\l 
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daida2{SXi ■ [X[ x [X21 x X'^] 



1_H^[2i|/Al |X2i|/Ai 



L21| 



da 



dz- 



2z 



-6X ■ [X' X [X' X X"]] = -2 In 



daSX ■ X".(6.2) 



Eq. (j3.12p should be kept in mind when deriving the above chain of calculations. As far as AS'ofr are concerned, 
one can see that the terms in Eq. (|5.8p arising upon contracting [Xi x X^]i[X2 x Xl^lj with X2iiX2ij are just the 
higher derivative terms which are not enhanced as InA^/^. They are displayed in Eq. (|5.10p . Using Eq. (|3.18p and 
(|5.9p for the short distance regularization and performing the manipulations analogous to the above, one can obtain 
the variations of the terms containing velocity X. Collecting all variations together and keeping only the terms 
corresponding to the lowest derivative ones in Eq. (|5.10p one gets 



SScS = - 



f $0 



d 



A. 



V47rAL 

which results in the equations of motion 

[XxX'] = 
Let us denote for the sake of brevity 



dtdaSX • <{ — [X' — [X X X']]ln-^ 
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2m 



X"ln- 



7 ^ 



X"ln^ 



[XxX'] X X']ln- 



dtda{6X[X X X']),(6.3) 



(6.4) 



h 

2m 



1 



and consider the large amplitude motions, when the nonlinearity in Eq. 
approximate solution obtained upon neglecting the terms cx l/c^ is 

[XxX'p =7Kn, XW =7Kb + l^'"^X' 



(6.5) 



could be essential. The zeroth order 



(6.6) 



where the first line in Eq. (|3.13p and the relation X' = [n x 6] are used. The longitudinal component of the velocity 



V^^'^'' cannot be determined from the first relation in Eq. (|6.6p . However, V^''"'' is locally unobservable, because of 
physical homogeneity of the vortex in the (local) longitudinal direction. It can be excluded by additional gauge choice 
XX' = 0. So, the zeroth order intrinsic velocity of the curved gauge vortex looks as 



.(0) 



X(o) 



bA|x' 

2m' 



In- 



(6.7) 



It is proportional to the curvature of the contour and is directed along the bi-normal vector. Its magnitude relative 
to the sound velocity Cs is 



|X(")| 



|X"|ln ^ 

2mCs c, 



^ln^«l, 
R £, 



where R = |X"|^^ stands for the radius of curvature of the contour. The large scale nonlinear motion of the contour 
is rather slow. One can see that the straight vortex does not possess an intrinsic motion as it should. 

To find the correction 0{l/c1) to the zeroth order solution it is necessary to obtain the time derivatives of the 
curvature k and torsion r as well as the triple of basic vectors of normal n, bi-normal b, and the tangent X'). Using 
Eq. (|3.13p . (|6.6p . the definitions k = (nX"), r = —{b'n) one obtains the equations 



dX' 
dn 



^{n'b — KTu), 
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and 



— = -7(2K'r + KT'), 



(6.9) 



Equations (j6.9p (however, without the factor 7) are known in hydrodynamics as Da Rios equations [l6j . Making use 
of Eq. (|3.13p . (|6.4p . (|6.6p . (|6.8p . and (|6.9p one can obtain the relation for finding the correction X*-^-* due to the term 
0(l/c?): 



[XWxX'l = -^In^ fn(K"-KT2) + 6(2K'r + KT')l 



(6.10) 



Using the above relation one gets 



^ ^ In 4^ r6(«:" - nr^) - n(2«'T + «;t')1 = jX^") In ^ In 4^ x 



217101 "' C 

[6(K"-KT2)-n(2K'T + Kr')] - 

The ratio of correction to the zeroth order velocity 



(6.11) 



As , Ar 



< 1 



is small but non-negligible, because the truly small ratio of the coherence length to the curvature radius squared is 
enhanced by the product of two large logarithms. 

Let us consider the small amplitude oscillations around the locally straight vortex. Then, locally, one may choose 
X' — 62, X{a^ t) = u{a, t) + e^fj, where u is orthogonal to the contour direction e^: • M(cr, t) = 0. Introducing the 
characteristic velocity 



X 



ln(AL/0 



and restricting oneself in Eq. 



ln(A,/0 ' 

by the linear approximation, one obtains the equation 
il ,, 1 . . , 

-rr — U H M X = 

Co 7 



(6.12) 



(6.13) 



for the two-dimensional vector u = {ux,Uy). The plain wave solution u^^y oc e*''^'^ '^^^ corresponds to the oscillation 
spectrum 



fc2 



and the relations among the amplitudes Uy = ±iux- In our case, j k /cq = 2k^ ln(AL/0 ln(As/0 ^ because the 
typical wave numbers k <^ 1/^, so that the square root can be expanded. Then the general solution of Eq. (j6.13p . 
assuming as usual the periodic boundary conditions, can be represented in the form of the superposition of the left 
and right circularly polarized waves: 



u{a, t) 



V2 



J2 (a/e-'^i'-HAe-'^^Oe^^ + cc. 



(6.14) 



I——0C 



Here, ki — 2ttI/L, L stands for the periodicity length, ^ = 0, 1, 2.. 



0J2 



fik} 
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2- "-Lin:^ 



2m C 
2 _ 2nog 



7 ' ?iln(A,/0 ' 2m ^ ' 



In—, 



(6.15) 
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and ai,3i being the arbitrary complex constants. The first branch ui coincides with the Friedel-De Gennes-Matricon 
waves [l3, [3 ' if O'^^ allows for the fact that here the scalar field mimics the condensate of the Cooper pairs with 
the mass m = 2me and charge q = 2e. It originates from the fact that the vortex segment being accelerated by the 
tension force oc X" acquires velocity Av cx n in the normal vector direction. The Magnus-like force oc [e^ x u] results 
in the circular motion of the segment in the plane perpendicular to and n, that is in the direction of bi-normal 
vector b = [X' x n]. The high frequency branch uj2 possesses the gap 

_ 2nog 
-Mn(A./0- 

It can be interpreted as follows. Since ojmin — ^2{ki = 0), this mode corresponds to the neglect of the term cx X" 
in the vortex equations of motion Eq. (|6.13p . Using Eq. (|5.12l) for the effective vortex mass together with Eq. (|4.8p . 
(|6.5p , and (j6.12p one can rewrite the equations of motion in this limiting case as 

^^^it = 2TTfino[ez x u], (6.16) 

which is the dissipation- and pinning- free part of the vortex equation of motion including the transverse force Q • 

The characteristic velocity cq in the case of gauge vortex exceeds the sound velocity Cs [see Eq. (|6.12p ] , hence both 
branches of small transverse oscillations are expected to be damped due to the emission of the Bogolyubov- Anderson 
excitations (sound waves in the large wave length limit). The calculation of the damping rate goes beyond the scope 
of the present paper. 



VII. THE CASE OF THE GLOBAL VORTEX. 



The global vortex can be treated in the same manner. The corresponding field theoretic model can be obtained 
from AHM one upon switching off the gauge field degrees of freedom and is called the Goldstone model [19] in particle 
physics and the Gross-Pitaevskii one iJj, i2i] in condensed matter physics. Omitting the chain of derivations completely 
analogous to those presented in proceeding sections let us write the effective action S'off = S'bg + AS't of the global 
model as the sum of the action of the background field S'bg and that of the integrated out fluctuations of the modulus 
and the phase of the scalar field ASf. Here Sh.y looks like 



hno 



ao - dtXh - -^{a + Vxb) 
Im 



(7.1) 



and Xb stands for the smooth background phase representing the possible potential fiow whose velocity is t^b = ^ Vxb, 
while ao — —dtXs and a = Vxs represent the time and space derivatives of the singular phase Xs responsible for the 
vortex. Their space Fourier components are given by the equations analogous to Eq. (|2.8p and (|2.7p . respectively: 



aok{t) = ~^^'^<^ / daak[Xa x X'Je 

a 



(7.2) 



From now on let us restrict ourselves by the single vortex with the unit quantum of circulation ria = 1. Then the 
coordinate space expressions found from Eq. (j7.2p by the inverse Fourier transform are 



ao(t, a;) 
a(i, x) 



da 
da 



{x~X)-[Xy. X'] 
\x-X\^ ' 
[X' X (a; - X)] 
\x-X\^ ' 



(7.3) 



X = X[t,a). The term originating from the integrated out fiuctuations can be obtained from Eq. (jS.ip by taking 
the limit of vanishing gauge coupling constant q and restoring the term oc a? /2m omitted when discussing the gauge 
vortex: 



h^riQ f d^k 



\Pk\ 



2m J (27r)3 ;i2^2/4 



(7.4) 
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where, for the sake of brevity, we introduce the quantity pk = J (Pxpe which is the Fourier amphtude of the 

quantity 



p = ao - dtXh 



2m 



{a + Vxb) 



(7.5) 



Notice that this p (taken at ao = 0, a = 0, that is, in the case of no vortex present) is just that combination of the 
derivatives of the Goldstone field Xb which was introduced in Ref. [13, [HI on the basis of the Gahlean invariance of 
the effective action. Since the typical momenta are subjected to the condition |fc| <C 1/^, where the healing length ^ is 
given by Eq. (|4.12p . the correction Eq. (|7.4p to the background action can be represented as the series in the number 
of derivatives 



A5f 



2mc1 



dtcfx 



1=0 



The leading order effective action in the case of no vortex present is, omitting the total derivative term, 

1 



2m 



r(9ab)'-(Vxb) 



(7.6) 



(7.7) 



which is just the action of the Goldstone mode in neutral superfluid. This agrees with the general analysis (2l|, where 
the coefficient in the effective action was fixed on the basis of the equations of motion of the Goldstone field, about 
which it is known that this field is the sound wave. 

Turning back to deriving the global vortex equations of motion let us set the smooth phase Xh to zero and neglect in 
Eq. (|7.6[) all terms (x . When finding the variation of the background effective action Eq. (|7.ip over X the following 
expressions derived from Eq. (17. 2p are necessary: 



(5ao(a;) 
5a{x) 



(2^ 
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{2nY 



dae-'^-^ <! [5X X X'] 



rx> 5X-[ky.X'].^ 

5X -[X -K X'] H 1_ liQ^ 



fc2 

k{k ■ [6X X X']) 

fc2 



(7.8) 



In fact, the time derivative term in Sa^ in Eq. (|7.8p vanishes when applied to the background action. The zeroth order 
equations of motion resulting from the variation of the background action are 



[X X X'] [a(X) X X'] = 0. 



(7.9) 



One can see from Eq. (|7.9p that upon neglecting the integrated out fluctuations of the phase and the modulus of the 
scalar field the classic law of the vortex motion 



-a(X) 



(7.10) 



is recovered, because ^a(X) is the velocity induced at the given vortex location x = X{a) by other elements of the 
vortex. Inserting here Eq. (j7.3p . using Eq. (j3.12p . and taking the low and upper integration limits to be, respectively, 
^ and R (see below), one obtains in the zeroth order 



X(o)( 



^1 



-a(X) 
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2m 



da- 



[X21 X X^] 



|^[XlxXnin| 



7s 



Kb, 



where b is the bi-normal vector, and 



n , R 



(7.11) 



(7.12) 



Recall that Xi 9 = X{ai^2), and k — \X"\. The result (|7.1ip coincides with the velocity of the vortex ring in 
superfluid Hell [1^. The large-scale motions of the contour in the case of the global vortex is slow, because 



|X(o)| 



hn R 
In — 

2mCs ^ 
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Here the curvature of the contour k is approximated by the inverse cutoff 1/R. 

As far as AS[ are concerned, we evaluate its contribution under the same assumption |fc| <^ 1/^ as made when 
discussing the local vortex case. However, in addition to the short distance divergence at r < ^ the energy per unit 
length of global vortex has the pertinent divergence at large distances due to the real gapless Goldstone field. As 
usual, one regularizes it by means of the integration cutoff at the distances r ^ R, where R is of the order of the 
curvature radius of the contour. We make here the large distance regularization by means of the large distance cutoff 
(fc^ + 1/R'^)~^, R —> oo. Performing the integration over momentum and making the same approximation as 
in the case of the local vortex one finds 

ASf « fdtda[[XxXfln^ + —{X-[X'xX"]A. (7.13) 

2mc- J [ C 4 J 

Contrary to the case of magnetic flux vortex, here the higher derivative term is multiplied by the large factor R^. The 
ratio of the second higher derivative term in the braces to the first lowest order derivative term is estimated with the 
help of Eq. (jTTTI) to be 

R^K^ 1 

< 1, 



lni?/C InR/C 

where we take the cutoff R ^ 1/k. So, the higher derivative here in the case of the global vortex is also suppressed, 
although without additional small factor Xj^K^ pertinent to the gauge vortex case. As in Sec. IVIl when obtaining the 
correction to the equations of motions, the higher derivative term is neglected. Analogously to the gauge vortex case, 
here the lowest derivative term in the action can be interpreted as the contribution of the kinetic energy of the vortex 
segment with the effective mass per length 

moff irh'^nQ R irfi^ R 

This coincides with the expressions obtained in Ref. 0, 0] , allowing for apparent typo (extra mass in the denominator) 
in Ref. [|. 

The correction to the zeroth approximation should be derived upon varying the nonlocal form of the action 
Eq. (j7.13p . Taking into account the zeroth order equations of motion, one finds 

[X X X'] = -[a{X) X X'] + -\\n^[dt[X x X'] x X']. (7.15) 

Recalling that a{X) w -lln-^b [see Eq. (|7.1ip ] Eq. (|7.15p looks similar to Eq. (|6.4p . Hence the correction due to 
excitations of the modulus and the phase of the scalar field can obtained from Eq. (|6.1ip by means of the replacements 
7 7g, \s,Xl R-- 

XH) ^ ^ In :^ [&(«;" _«;r2)-n(2«:'r + «;t')] = 
2mcj 4 ■' 

|X(o)|^2 ^^2 R _ ^^2) _ ^(2k't + kt')1 -. (7.16) 

4 K 

The ratio of the correction to the zeroth order approximation is small, 

— t, K m — < 1, 



|X(o)| e 

but non-negligible, because of the large logarithm squared. As in the case of the magnetic vortex, the large scale 
nonlinear motion is slow. 

The linear small amplitude motion of the global vortex is analyzed by taking the limit of locally straight contour, 
see preceding section for more detail. The corresponding equations of motion looks similar to Eq. (|6.13p . in which 
one should make the replacements cq ^ Cg, because the ratio of logarithms in Eq. (|6.12p drops in the global vortex 
limit, and J ^ "fg- Hence, the small transverse motions of the global vortex are characterized by the velocity of sound 
which means that the kinematics is exactly on the threshold of the emission of sound waves, and the phase space of 
the emitted waves shrinks to zero. By this reason the small amplitude motions are not expected to be damped. The 
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oscillation spectrum is given by Eq. (|6.15p with the proper replacements just mentioned, together with Al, As R. 
The lower branch of the spectrum is well-known as the Kelvin waves [ll, [T^, [13 . The branch with the gap 

_ 2nog 
h\n{R/0 

is interpreted similarly to that of the local case as being due to the acceleration of the vortex segment with the mass 
per length Eq. (|7.14p subjected to the Magnus force. The equations of motion look the same as Eq. (|6.16p . with the 
effective mass Eq. (|5.12p replaced by one in Eq. (|7.14[) . 



VIII. DISCUSSION. 



The usual attitude to the effective lagrangians is that they do not refer to underlying basic theory and its fun- 
damental constituents but are formulated in terms of effective degrees of freedom. In the case of condensed mater 
physics they are the modulus and phase of the order parameter together with the screened electromagnetic field in 
case of charged superfluid, and the same except the electromagnetic fields in case of neutral superfluids. Guided by 
analogous approach based on the field theoretic models possessing the property of spontaneous symmetry breaking 
we obtain here the effective actions and the equations of motion to the lowest number of derivatives over contour 
parameter and time for the local (gauge) and global vortices. Of course, the models used in the present paper are 
idealized in that they ignore the friction and pinning effects, but the inertial properties of the vortex state are most 
conveniently treated in the framework of such models. The key feature of the present consideration is the usage of 
the explicitly transverse electric gauge field strength in the study of the gauge vortex dynamics. If the field excita- 
tions were not included, the resulting effective action of the background field configuration would contain divergence 
at large distance. The crucial role of the excitations is that they cancel this divergence so that the gauge vortex 
dynamics becomes finite at large distances as it should, because all gauge fields in spontaneously broken gauge model 
are screened at the distances larger than the penetration depth. As a consequence, the higher derivative terms in 
the effective action for curved gauge vortex are small. Taking into account the scalar field excitations in the global 
models like the Gross-Pitaevskii or the Goldstone one results in the correct Galilean-invariant form of effective action 
for the Goldstone mode suggested in Ref. [2l[ on the symmetry grounds. 

I am grateful to G. I. Shvetsova for encouragement and warm support. 
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